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Abstract 

The behavior of Laplacian eigenfunctions in domains with branches is 
investigated. If an eigenvalue is below a threshold which is determined 
by the shape of the branch, the associated eigenfunction is proved to 
exponentially decay inside the branch. The decay rate is twice the square 
root of the difference between the threshold and the eigenvalue. The 
derived exponential estimate is applicable for arbitrary domains in any 
spatial dimension. Numerical simulations illustrate and further extend 
the theoretical estimate. 
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1 Introduction 

The Laplace operator eigenfunctions play a crucial role in different fields 
of physics: vibration modes of a thin membrane, standing waves in optical 
or acoustical cavity resonators, the natural spectral decomposition basis for 
diffusive processes, the eigenstates of a single trapped particle in quantum 
mechanics, etc. For the unit interval, the spectrum of the Laplace operator 
is particularly simple, and the eigenfunctions are just linear combinations 
of Fourier harmonics {e tlvnx }. The oscillating character of eigenfunctions 
is then often expected for domains in two and higher dimensions. The 
"scholar" examples of the explicit eigenbases in rectangular and circular 
domains strongly support this oversimplified but common view. At the 
same time, the geometrical structure of eigenfunctions may be extremely 
complicated even for simple domains (e.g., the structure of nodal lines of 
degenerate eigenfunctions in a square [1]). The more complex the domain 
is, the more sophisticated and sometimes unexpected the behavior of the 
associated eigenfunctions may be. For instance, numerous numerical and 
experimental studies of the eigenvalue problem in irregular or (pre)fractal 
domains revealed the existence of weakly localized eigenfunctions which have 
pronounced amplitudes only on small subregions of the domain [2-13]. 

In this paper, we investigate the behavior of Laplacian eigenfunctions in 
domains with branches. We show that certain eigenfunctions are "expelled" 
from the branch, i.e., their amplitude along the branch decays exponentially 
fast. A similar "expulsion" effect is well known in optics and acoustics: a 
wave of wavelength i cannot freely propagate inside a rectangular channel 
of width b smaller than 1/2 because of the exponential attenuation ~ e ~ xn / b 
along the channel [2, 14]. We extend this classical result to arbitrary domains 
with branches of arbitrary shape. We derive a rigorous exponential estimate 
for the Z/2-norm of the eigenfunction in cross-sections of the branch. We 
obtain the sharp decay rate which generalizes and refines the classical rate 
n/b. Although this problem is remotely related to localization of waves in 
optical or acoustical waveguides (e.g., infinite bended tubes [14-16]), we 
mainly focus on bounded domains. 

It is worth noting that the exponential decay of eigenfunctions of Schrodinger 
operators in free space (so-called strong localization) has been thoroughly 
investigated in physical and mathematical literature (see, e.g., [17-21]). The 
first exploration of this problem for arbitrary Schrodinger potential bounded 
from below was given by Schnol' [17] who proved an exponential decay of 
eigenfunctions in which the decay rate was related to the distance between 
the corresponding eigenvalue and the essential spectrum. This result is of 
remarkable generality because the essential spectrum may be arbitrary, for 
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example with gaps, and may not consist of positive axis. A sharp esti- 
mate for the decay rate was made by Maslov who reduced the problem to 
a differential inequality [20, 21]. Anderson discovered the exponential de- 
cay of eigenfunctions of the Schrodinger operator with random potentials 
[22]. This phenomenon, known as the Anderson localization, has been in- 
tensively investigated (see reviews [23, 24]). Although the Laplace operator 
in a bounded domain is a much simpler mathematical object, the properties 
of its eigenfunctions are still poorly understood. 

The paper is organized as follows. In Sec. 2, we start by considering a 
two-dimensional domain with a rectangular branch. In this special case, the 
estimates are derived in a rather elementary and straightforward way that 
helps to illustrate many properties of eigenfunctions. Section 3 presents the 
analysis for domains with branches of arbitrary shape in any spatial dimen- 
sion. We provide a sufficient condition on the eigenvalue, under which the 
related eigenfunction decays exponentially inside the branch. Sec. 4 presents 
numerical examples which illustrate the theoretical results and suggest new 
perspectives for further investigations. The paper ends by conclusions, in 
which the main results are summarized and their consequences are discussed. 

2 Rectangular branch 

We begin with the following example. We consider the Dirichlet eigenvalue 
problem 

(1) — An = An, (x,y) 6 D, u\qd = 

in a planar bounded domain D = V U Q which is decomposed into a basic 
domain V of arbitrary shape and a rectangular branch Q with sides a and 
b: Q = {(x,y) G M 2 : < x < a, < y < b}, as illustrated on Fig. la. 
We assume that the eigenvalue A is smaller than the first eigenvalue of the 
Laplace operator in the cross-section of the branch (i.e., the interval [0, b]): 

(2) A < ir 2 /b 2 . 

Under this condition, we aim to show the exponential decay of the associated 
eigenfunction n in the branch Q. 

In this illustrative example, the derivation is based on an explicit upper 
bound for the norm of the eigenfunction n. One can easily check that a 
general solution of Eq. (1) in the rectangular branch Q has a form 

oo 

(3) n(x, y) = ^2 c n sinh(7„(o - x)) sm(irny/b), 

n=l 
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Figure 1. A bounded domain D is the union of a basic 
domain V of arbitrary shape and a branch Q. Our goal 
is to estimate how an eigenfunction of the Laplace opera- 
tor in D decays inside the branch, (a) rectangular branch 
Q = {(x,y) £l 2 : < x < a, < y < b} (the energetic 
norm is estimated in the subdomain Q(xq) shown by shad- 
owed region); (b) branch of arbitrary shape [the Z^-norm is 
estimated in the cross-section fi(x)]; (c) three-dimensional 
increasing branch. 



where j n = (f-n) 2 — A, and c n are constants. We consider the energetic 
norm of the function u in the subdomain Q(xq) = {(x, y) E M 2 : xq < x < 



a, < y < b} which is defined as 

( 4 ) W Vu \\l 2 (Q{x )) 



(Vu, Vu)dxdy. 



Substituting Eq. (3) to (4) yields 



|Vu| 



L 2 (Q(x )) 



(5) 



n=l 



, 7T 



{'^n) 2 smh 2 {-f n (a - x)) 
b 



■Co 



+ 7^ cosh 2 (7 n (a - x)) 



dx. 



Using elementary inequalities for the integral (see Appendix A), one gets 

oo 

(6) HVullLwCxo)) ^ Ce^^nc 2 sinh 2 ( 7n a), 



n=l 
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where C is an explicit constant. The trace theorem provides the upper 
bound for the above series (see Appendix A) 

\\^\\l (Qixo)) <Cie^, 

with another explicit constant C\. So, we established the exponential decay 
of the energy 1 1 Vu\ | 2 in the branch Q with the decay rate 271 = I^J (j) 2 — A. 

From this estimate it is easy to deduce a similar estimate in the L2-norm, 

(7) IMlLw-o)) 25 / u 2 dxdy<C 2 e- 2 ^, 

where C2 is another constant. The derivation implies that the estimate is 
sharp, i.e. the decay rate cannot be improved in general. 

It is worth noting that no information about the basic domain V was used. 
In particular, the Dirichlet boundary condition on dV can be replaced by 
arbitrary boundary condition under which the Laplace operator is still self- 
adjoint. 



3 Branch of arbitrary shape 

Now we show that the above estimate remains valid for eigenfunctions in 
a much more general case with a branch of arbitrary shape in R n+1 (n = 
1, 2, 3, ...). We consider the eigenvalue problem 

(8) - Au(x, y) = Xu(x, y) (i,y)efl, u\ dD = 0, 

where D = V U Q is decomposed into a basic bounded domain V C M n+1 
of arbitrary shape and a branch Q C R ra+1 of a variable cross-section profile 
Q(x) C R n (Fig. lb,c): 

Q = {( x ,y) G R n+1 : y G Q(x), < x < a}. 

Each cross-section Q(x) is a bounded domain which is parameterized by x 
from to a. The boundary of the branch Q is assumed to be piecewise 
smooth [25]. Although weaker conditions on the boundary could potentially 
be used, their justification would require a substantial technical analysis 
which is beyond the scope of the paper. 

For a fixed x, let (J>i(x) be the first eigenvalue of the problem 

(9) -A x <Ky) = MiW(y) yefi(4 <fe = o, 

where Aj_ is the n-dimensional Laplace operator. We denote 

(10) n= inf m(x) 

0<x<a 
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the smallest first eigenvalue among all cross-sections of the branch. For 
example, if Q(x) = [0, b] (independent of x), one has /x = n 2 jb 2 and retrieves 
the example from Sec. 2. 

Now, we can formulate the main result of the paper. If the basic domain 
V is large enough so that 

(11) A < fi, 

we prove that the squared L2-norm of the eigenfunction u in the cross-section 



(12) I(x )= J u 2 (x ,y)dy, 

Q{xq) 

decays exponentially with xq: 

(13) I(x ) < /(0)e- /3x ° (0 < x < a), 
with the decay rate 



(14) = V2 y/fi-X. 

Moreover, if the branch profile Q(x) satisfies the condition 

(15) (ex,n(x,y))>0 V(x,y)G9Q, 

where e x is the unit vector along the x coordinate, and n(x,y) the unit 
normal vector at the boundary point (x, y) directed outwards the domain, 
then the decay rate is improved: 



(16) P = 2vV-A. 

Qualitatively, the condition (15) means that the branch Q is not increasing. 
Note that the example of a rectangular branch from Sec. 2 shows that the 
decay rate in Eq. (16) is sharp (cannot be improved). 

The proof consists in three steps. 

(i) First, we derive the inequality 

(17) I"(x ) > c 7 2 I(x ), 

where c = 2 for arbitrary branch and c = 4 for branches satisfying the con- 
dition (15), and 7 = y/jl — A. This type of inequalities was first established 
by Maslov for Schrddinger operators in free space [20, 21]. For this purpose, 
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we consider the first two derivatives of I(xq): 

du 
dx 



(18) I'(x ) = 2 / u^dy, 



n(x ) 

d 2 u , „ f ( du\ 2 



(19) />„) . + 

n(a; ) H(a;o) 

where the boundary condition u^g = cancels the integrals over the "lat- 
eral" boundary of Q{xq). 

The first integral in Eq. (19) can be estimated as 
f d 2 u f 

/ u ~Q^i dy = / u [Au-A ± u]dy 

(20) n ^ Q ^ 

= J (V_i_«, V ± u)dx -A y u 2 dy > (ji — X) J u 2 dy, 

£1(xq) £1{%o) tt(xo) 

where we used the Friedrichs-Poincare inequality for the section fJ(xo) ( see 
Appendix B.l): 

(21) J (V ± u,V±u)dy >in(x ) J u 2 dy>fi J u 2 dy, 

and Hi(xo) > [/, by definition of \i in Eq. (10). Since the second term in Eq. 
(19) is always positive, one has 

I"{x )>2 J u-Tj^dy > 2(/i - A) J u 2 dy, 
Q(x ) ft(x ) 
from which follows the inequality (17) with c = 2. 

If the condition (15) is satisfied, a more accurate estimate of the second 
term in Eq. (19) follows from the Rellich's identity (see Appendix B.2): 

2 



/ (£) dy= J W ± u,v ± u)dy-x j u2d y 

n(x ) ci(x ) a(x ) 

r (n \ 2 

+ J (£j (exMS))dS, 



8Q(x )\fl(x ) 

where Q(xq) denotes the "right" part of the branch Q delimited by £1(xq): 
(23) Q(x ) = {(x,y) G R n+1 : y G Q(x), x < x < a} . 
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The condition (15) implies the positivity of the last term in Eq. (22), which 
can therefore be dropped off in order to get the following estimate: 

/ (£) dy ~ / ( V ^' V ^ n )^- A / ^ 2 dy>(u-X) J u 2 dy. 

Combining this result with (20), one gets the inequality (17) with c = 4. 

(ii) Second, we prove the following relations 

(24) 1(a) =0, l'(a) = 0, I(x ) ± 0, l'(x ) < 

for < xq < a. In fact, taking into account Eq. (18) and applying the 
Green's formula in subdomain Q(xq), one obtains 

(25) l'(x ) = -2 J uAudxdy-2 J (Vu,Vu)dxdy, 

Q(x ) Q(x ) 

where the boundary condition u\qq^ Xo ^\q^ X( ^ = canceled boundary inte- 
grals. The second term can be estimated by using again the Friedrichs- 
Poincare inequality (21): 

2 



(26) 



(Vu, Vu)dxdy = 

Q(x ) Q(x ) 



du\' 



dxdy 



a a 

> j dx J (y±u, V ±u)dy > \x J dx J u 2 dxdy = [i J u 2 dxdy. 



xo fl{x) x Q(x) Q(xo) 

The equation Au = —Xu yields then 



-I'(x ) = -2X J u 2 dxdy + 2 J {Vu,Vu)dxdy 

Q(xo) Q(xo) 

> 2(u - X) J u 2 dxdy > 0, 



Q(xo) 

so that I(xq) monotonously decays. We emphasize that there is no request 
on a monotonous decrease of any kind for the cross-section O(x). 

Sending xo to a in Eq. (25), one gets I'(a) = 0. 

Finally, we show that I(xq) ^ for < xo < a. Indeed, if there would 
exist xo such that I(xq) = 0, then the restriction of u to the subdomain 
Q(xq) is a solution of the eigenvalue problem in Q(xq): 

-Au = Xu (x,y)£Q(x ), u\dQ( X0 ) = °- 
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Multiplying this equation by u and integrating over the domain Q(xq) lead 
to 



A j u 2 dxdy = J (Vtt, Vu)dxdy. 



Q(x ) Q(x ) 

On the other hand, the Friedrichs-Poincare inequality (26) yields 

J (Vu, Vujdxdy > fi J u 2 dxdy, 
Q(x ) Q{xo) 
from which A > fi, in contradiction to the condition (11). 

(iii) Third, we establish the exponential decay (13) following the Maslov's 
method [20, 21]. The multiplication of the inequality (17) by I'(xq) yields 

((J') 2 )' < c 7 2 (J 2 )'. 

After integrating from xq to a, one gets 

-(/'(*o)) 2 < -ci 2 I 2 (x ), 

where 1(a) = I 1 (a) = from Eq. (24). Taking in account that /' < 0, one 
deduces 

-I'(x ) > V^jl(x ). 
Dividing by I(xq) and integrating from to xq lead to Eq. (13) that com- 
pletes our proof. 

As in Sec. 2, no information about the basic domain V was used so 
that the Dirichlet boundary condition on dV can be replaced by arbitrary 
boundary condition under which the Laplace operator in D remains self- 
adjoint. 



4 Discussion 



The main result of Sec. 3 states that if an eigenvalue A is below the smallest 
eigenvalue [i in all cross-sections f2(x), the associated eigenfunction u decays 
exponentially in the branch Q. This condition on the eigenvalue A in D can 
be replaced by another condition 

(27) k < fi 
on the eigenvalue k in the basic domain V: 

(28) -Acf) = of, (x, y) G V, <P\ dv = 0. 

According to the Rayleigh's principle, the inequality (27) is sufficient for the 
existence of the eigenvalue A lying below \x (see Appendix B.l). Although 
the condition (27) is stronger than (11), it may be preferred for studying 
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n 


2a 


2c 


2d 


2e 


2f 


2g 


2h 


2i 


2j 


1 


19.33 


19.66 


19.39 


19.64 


19.58 


19.39 


19.33 


19.39 


19.39 


2 


47.53 


48.95 


47.58 


48.94 


48.59 


47.86 


47.54 


47.58 


47.58 


3 


49.32 


49.35 


49.33 


49.35 


49.33 


49.32 


49.32 


49.33 


49.33 


4 


78.83 


78.75 


77.85 


78.90 


78.54 


78.85 


78.84 


77.85 


77.85 


5 


93.12 


97.87 


94.41 


97.69 


97.09 


94.48 


93.12 


94.40 


94.41 


6 


98.70 


98.71 


98.67 


98.71 


98.66 


98.71 


98.71 


98.67 


98.67 


7 


126.1 


127.8 


125.2 


127.9 


127.3 


126.9 


126.1 


125.2 


125.2 


8 


128.0 


128.3 


128.0 


128.3 


128.0 


128.1 


128.0 


128.0 


127.7 


9 


151.7 


166.1 


154.1 


165.8 


164.6 


158.5 


151.6 


150.4 


128.0 


10 


167.7 


167.8 


167.8 


167.8 


167.8 


167.7 


167.7 


156.8 


153.6 


11 


167.8 


177.5 


176.5 


177.1 


177.0 


175.8 


171.8 


167.8 


167.8 


12 


175.3 


193.9 


182.3 


196.9 


183.4 


196.5 


176.1 


176.7 


167.8 


13 


191.7 


196.1 


196.1 


197.5 


195.3 


197.4 


196.4 


185.3 


176.7 


14 


196.4 


197.5 


197.4 


239.6 


197.4 


229.6 


197.4 


197.1 


185.4 


15 


197.4 


245.1 


229.1 


244.4 


244.1 


245.9 


204.5 


197.4 


195.2 


16 


218.8 


246.8 


246.0 


246.8 


246.2 


250.7 


235.5 


218.5 


195.5 


17 


245.7 


254.3 


249.2 


254.1 


252.4 


256.7 


245.8 


242.8 


197.4 


18 


246.7 


256.7 


256.6 


256.7 


256.6 


284.6 


250.5 


246.1 


214.8 


19 


252.5 


284.7 


269.6 


285.8 


259.3 


285.5 


256.7 


250.6 


229.0 


20 


256.6 


286.3 


285.9 


286.3 


283.2 


304.1 


278.7 


256.6 


245.9 



Table 1. First 20 eigenvalues of the Laplace operator in do- 
mains shown on Fig. 2. 



different branches Q attached to the same domain V (in which case, the 
eigenvalue k has to be computed only once). 

Looking at the derivation of the inequality (13), many questions naturally 
appear: How accurate the exponential estimate is? Is the condition (15) 
necessary for getting the sharp decay rate according to Eq. (16)? How do 
the eigenfunctions with A larger than fx behave? Is the exponential decay 
applicable for other boundary conditions? In the next subsection, we address 
these questions through numerical simulations in planar domains. 

4.1 Two-dimensional branches 

In order to answer these and some other questions, we consider several planar 
bounded domains D which are all composed of the unit square V as the 
basic domain and a branch Q of different shapes (Fig. 2). The eigenvalue 
problem (8) in these domains with Dirichlet boundary condition was solved 
numerically by Matlab PDEtools. Once the eigenfunctions and eigenvalues 



A. DELITSYN, B.-T. NGUYEN, D. S. GREBENKOV 

















a 



(a) 1 



•I'D 



(b) 





x 0< ; 1 

)b/2 




a 



x 



]6/2 



(d) 




Figure 2. The unit square (basic domain V) and several 
shapes of the branch Q: (a) rectangular branch; (b) half 
of the domain 'a'; (c) narrow-than-wide channel; (d) wide- 
than-narrow channel; (e) increasing branch with the width 
linearly changing from 6/2 to b; (f) branch with a partial 
cut at the middle; (g) tilted rectangular branch; (h) circular 
branch; (i) branch with a small broadening in the middle; 
(j) bifurcating branch. We set a = 1 and b = 1/4 in all 
cases, except 'g' and 'h', for which a = l/v2 and a = 5/8, 
respectively. For shapes 'c', 'd', V, 'f, 'i' and 'j', the branch 
is up shifted by 1/8 in order to break the reflection symmetry 
(for cases 'g' and 'h', there is no shift because the branch itself 
has no reflection symmetry). 
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are found, we approximate the squared L2-norm of the eigenfunction u n in 
the subregion Q(xq) = {(x,y) G Q : xq < x < a} of the branch Q as 

(29) J n (x )= f ul(x,y)dxdy ^y2^p-^2u 2 n (xj,yj), 

Q(xo) 1 J ~ L 

where the sum runs over all triangles T of the mesh, S(T) being the area of 
the triangle T, and {(xj, yj)}j=i,2,3 its three vertices. Since the eigenfunc- 
tions are analytic inside the domain, the error of the above approximation is 
mainly determined by the areas of triangles. In order to check whether the 
results are accurate or not, we computed the function J n (xo) at different 
levels k of mesh refinement (once the initial triangular mesh is generated 
by Matlab, each level of refinement consists in dividing each triangle of the 
mesh into four triangles of the same shape). Figure 3 shows the resulting 
curves for the rectangular branch (Fig. 2a). For the first eigenfunction, all 
the curves fall onto each other, i.e. Ji(xo) is independent of k, as it should 
be. In turn, the curves for n = 3 coincide only for small xq but deviate from 
each other for larger xq. The higher k, the closer the curve to the expected 
exponential decay. This means that even 6 levels of mesh refinement (i.e., a 
mesh with 438272 triangles) is not enough for an accurate computation of the 
integral ^(xo). Among the 20 first eigenfunctions, similar deviations were 
observed for n = 3,4, 8, 10, 14, 15, 17. The specific behavior of these eigen- 
functions seems to be related to their reflection symmetry. Although several 
improvements could be performed [e.g., higher-order integration schemes in- 
stead of Eq. (29)], this is too technical and beyond the scope of the paper 
because our numerical computations aim only at illustrating the theoretical 
estimates. After all, the deviations become distinguishable only in the re- 
gion of xo for which J n (x$) is negligible. For all data sets discussed below, 
we checked the accuracy by performing computations with different k and 
presented only the reliable data with k = 5 (such meshes contain between 
100000 and 170000 triangles, except for the case on Fig. 2f with 671744 
triangles). 

The exponential decay (13) for I(xo) implies that of J n (xo): 

a a 

(30) J n (x ) = J dx I n (x) < J dx l n (0)e~ 2 ^ x < I01 e - 2 ^ x \ 

XO XO 

where 

and we take c = 4 even if the sufficient condition (15) is not satisfied. In 
what follows, we will check numerically the stronger inequality 

(31) Jn(x ) < J n (0)e- 2 ^ X0 , 
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10' 




10 







0.2 



0.4 



0.6 



0.8 



Figure 3. (Color online) Computation of J n (xo) at different 
levels k of mesh refinement of the rectangular branch on Fig. 
2a: k = 3 (symbols, 6848 triangles in the mesh), k = 4 (solid 
lines, 27392 triangles), k = 5 (dashed lines, 109568 triangles) 
and k = 6 (dash-dotted lines, 438272 triangles). For n = 1, 
all these curves fall onto each other, confirming the accurate 
computation which is independent of the mesh size. In turn, 
the curves for n = 3 coincide only for small xq but deviate 
from each other for larger xq. The higher k, the closer the 
curve to the expected exponential decay (Fig. 4). 

from which (30) follows, because 



It is worth stressing that the inequality (31) for the squared L2-norm 
J n (xo) in the subregion Q(xq) is a weaker result than the inequality (13) 
for the squared L2-norm I n (xo) in the cross-section CI(xq). However, the 
analysis of I n (xo) would require an accurate computation of the projection 
of an eigenfunction u n , which was computed on a triangular mesh in D, onto 
the cross-section Q(xa). Since the resulting I n (xo) would be less accurate 
than J ti {xq), we focus on the latter quantity. 

Rectangular branch 

We start with a rectangular branch of width 6 = 1/4 (Fig. 2a), for which 
fJL = TT /(1/4) 2 ~ 157.91..., and there are 9 eigenvalues A n below fj, (the 
first 20 eigenvalues are listed in Table 1). Figure 4a shows J n (xo) for four 
eigenfunctions with n = 1,8,9,13. These four eigenfunctions are chosen 



a 



a 
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x o 




Figure 4. (Color online) The squared L2-norm, J n (xo), of 
four eigenfunctions with n = 1,8,9,13 (symbols) for the 
rectangular branch on Fig. 2a. The estimate (31) with 
H = 7r 2 /(l/4) 2 is plotted by dashed (n = 1), solid (n = 8) and 
dash-dotted (n = 9) lines. The estimate with /i = 7r 2 /(l/8) 2 
is shown for n = 8 by dotted line. 

to illustrate different possibilities. For the eigenmodes with n = 1,2,5,7,9 
(illustrated by n = 1,9), the estimate is very accurate. In this case, the 
decay rate 2^f n is sharp and cannot be improved. There is another group 
of eigenfunctions with n = 3, 4, 6, 8 (illustrated by n = 8), for which J n {xo) 
is significantly smaller than the estimate. For xo < 0.4, J${xq) decays as 
J§(0) exp[— 27grr ], where 2^ = 2-^4/i — As is the improved decay rate (for 
larger xq, the computation is inaccurate as explained earlier and its result 
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is not shown). In order to understand this behavior, one can inspect the 
shape of the eigenfunction us(x,y) shown on Fig. 4b. This function is anti- 
symmetric with respect to the horizontal line which splits the domain D 
into two symmetric subdomains. As a consequence, us(x,y) is along this 
line and it is thus a solution of the Dirichlet eigenvalue problem for each 
subdomain (Fig. 2b). The width of the branch in each subdomain is twice 
smaller so that one can apply the general estimate with // = 4/i. This is a 
special feature of all symmetric domains. 

If the branch was shifted upwards or downwards, the reflection symmetry 
would be broken, and the decay rate 27^ would not be applicable any more. 
However, if the shift is small, one may still expect a faster exponential 
decrease with the decay rate between 27 n and 27^. This example shows that 
the estimate (13) may not be sharp for certain eigenfunctions. At the same 
time, it cannot be improved in general, as illustrated by the eigenfunctions 
with n = 1,9. 

The last curve shown on Fig. 4a by pluses, corresponds to n = 13, for 
which A13 > /i, and the exponentially decaying estimate is not applicable. 
One can see that the function Jisixo) slowly varies along the branch. This 
behavior is also expected from its shape shown on Fig. 4b. 

Narrow/wide and wide/narrow branches 

In two dimensions, any cross-section £l(x) of a branch is a union of intervals. 
If £(x) is the length of the largest interval in f2(x) then the first eigenvalue 
fi(x) in Q(x) is simply 7r 2 /£(x) 2 . Whatever the shape of the branch is, the 
bound fi for the exponential decay is then set by the length of the largest 
cross-section b = max l(x) according to [i = it 2 /b 2 . At first thought, this 

0<x<a 

statement can sound counter-intuitive because one could expect that the 
asymptotic behavior would be determined by the smallest cross-section. In 
order to clarify this point, we consider an example on Fig. 2c for which 
fl = 7T /(1/4) 2 . Although the narrow channel of width 6/2 strongly attenu- 
ates the amplitude of the eigenfunction, it does not imply the exponential 
decay with a hypothetical rate 2y/4fi — A along the next wider channel (Fig. 
5). For instance, the function Ji(xo) exponentially decays with the rate 
2-v/4/i — A (dashed line) up to xq 0.45. However, this decay is slowed 
down for x > 0.45. The theoretical estimate with the decay rate 2y//j, — A 
(solid line) is of course applicable for all xo but it is not sharp. 

If the wide channel is placed first (as shown on Fig. 2d), the situation is 
different. The decay rate 2-y/ /i — A, which is set by the largest cross-section 
length b = 1/4, can be used in the wide part of the branch. Although 
this result is also applicable in the narrow part, the estimate here can be 
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Figure 5. (Color online) The squared L2-norm, J u {xq), of 
three eigenfunctions with n = 1,9,12 (symbols) for the 
narrow-than-wide branch on Fig. 2c. The estimate (31) with 
fx = 7r 2 /(l/4) 2 is plotted for n = 1 by solid line. The hypo- 
thetical estimate with fi = it 2 / (1/8) 2 is shown by dash-dotted 
(n = 1) and dashed (n = 9) lines. The vertical dotted line 
indicates the connection between two parts of the branch. 



improved. In fact, if one considers the unit square and the wide branch as 
a basic domain, the decay rate for the narrow branch, 2\/4/x — A, is the set 
by its width 6/2 = 1/8 (Fig. 6). More generally, if the branch is a union of 
several branches with progressively decreasing widths, one can combine the 
exponential estimates with progressively increasing decay rates. 
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Figure 6. (Color online) The squared L2-norm, J n (xo), of 
three eigenfunctions with n = 1, 3, 12 (symbols) for the wide- 
than-narrow branch on Fig. 2d. The estimate (31) with 
fx = 7r 2 /(l/4) 2 is plotted by dotted lines for n = 1 and n = 3. 
The combined estimate (for the wide part with n = 7r 2 /(l/4) 2 
and for the narrow part with \i = 7r 2 /(l/8) 2 ) is shown by 
solid (re = 1), dashed (re = 3) and dash-dotted (n = 12) lines. 
The vertical dotted line indicates the connection between two 
parts of the branch. 



Increasing branch 



In the previous example on Fig. 2c, the supplementary condition (15) was 
not satisfied on a part of the branch boundary. Nevertheless, the numerical 
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Figure 7. The squared L2- n orm, J u (xq), of two eigenfunc- 
tions with n = 1,8 (symbols) for the increasing branch on 
Fig. 2e. The estimate (31) with fj, = 7r 2 /(l/4) 2 is plotted by 
solid (n = 1) and dash-dotted (n = 8) lines. The hypotheti- 
cal estimate with \i = 7r 2 /(l/8) 2 is shown by dashed (n = 1) 
and dotted (n = 8) lines. 

computation confirmed the sharp decay rate (16). In order to check the 
relevance of the condition (15), we consider a linearly increasing branch 
shown on Fig. 2e. Although the condition (15) is not satisfied at any point, 
the sharp decay rate (16) is again applicable, as shown on Fig. 7. In a future 
work, one may try to relax the condition (15) or to provide counter-examples 
showing its relevance. 

Branch with a cut 

A small variation in the shape of the domain is known to result in small 
changes in eigenvalues of the Laplace operator, at least at the beginning of 
the spectrum (see Table 1). In turn, the eigenfunctions may be very sensitive 
to any perturbation of the domain. Bearing this sensitivity in mind, one can 
check the robustness of the exponential decay. We consider a horizontal cut 
at the middle of the rectangular branch (Fig. 2f). If the cut went along the 
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whole branch, it would be equivalent to two separate rectangular branches of 
width 6/2. In this case, the theoretical estimate could be applied individually 
to each branch, and the value \x = 7r 2 /(6/2) 2 would be 4 times larger than 
for the original rectangular branch of width b. For a partial cut, whatever its 
length is, the theoretical decay rate is again determined by \x = 7r 2 /(l/4) 2 
as for the rectangular branch. At the same time, it is clear that the cut 
results in a stronger "expelling" of eigenfunctions from the branch. 

Figure 8a shows the squared L2-norm, J n (xo), of three eigenfunctions 
with n = 1, 11, 12. For the first eigenfunction, we plot two exponential 
estimates, one with the rigorous value fi = 7r 2 /(l/4) 2 and the other with 
a hypothetical value // = 7r 2 /(l/8) 2 = 4/i for a twice narrower branch (if 
the cut was complete). Although the first estimate is of course applicable 
for the whole region, it is not sharp. In turn, the second estimate is sharp 
but it works only up to xq < 0.4. One can see that for xo > 0.4, the slope 
of In Ji(xo) is given by the first estimate. The best estimate would be a 
combination of these two but its construction depends on the specific shape 
of the branch. It is worth noting that the transition between two estimates 
(i.e., the point 0.4) does not appear at the end of the cut (here, at 0.5). 
This means that one cannot consider two subregions (with and without cut) 
separately in order to develop the individual estimates. 

The 11th eigenvalue An « 177.0 exceeds [i so that the rigorous decay 
rate is not applicable. In turn, the use of the hypothetical value ji 1 provides 
a good estimate up to xq < 0.3 but fails for xq > 0.3. Finally, the 12th 
eigenfunction has no exponential decay (in fact, it is localized in the branch). 

Tilted and circular branches 

Another interesting question is the parameterization of the branch Q. In 
Sec. 3, a variable shape of the branch was implemented through the cross- 
section profile n(x) where x ranged from and a. The choice of the x coordi- 
nate is conventional and any other coordinate axis could be used instead of x 
(by rotating the domain). However, the freedom of rotation may lead to inac- 
curate estimate for the decay rate. This point is illustrated by the domain on 
Fig. 2g with a rectangular branch tilted by 45°. Applying formally the esti- 
mate (13), one expects the exponential decay with /i = 7r 2 /(l/4) 2 « 157.91. 
At the same time, if the whole domain was turned clockwise by 45° [or, 
equivalently, if the branch was parameterized along the axis x' in the di- 
rection (1, 1)], the decay rate would be set by \i = it 2 /(1/4/V2) 2 » 315.83. 
It is clear that the behavior of eigenfunctions does not depend neither on 
rotations of the domain, nor on the parameterization of the branch. Figure 
9 confirms the exponential decay with the latter threshold \x. 
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Figure 8. (Color online) The squared L2-norm, J n (xo), of 
three eigenfunctions with n = 1, 11, 12 (symbols) for the 
branch with a cut on Fig. 2f. The estimate (31) with 
fj, = 7r 2 /(l/4) 2 is plotted by solid (n = 1) and dash- 
dotted (n = 11) lines. The hypothetical estimate with 
fj, = 7r 2 /(l/8) 2 for n = 1 is shown by dashed line. The 
vertical dotted line indicates the end of the cut. 

A circular branch on Fig. 2h is another example, for which the choice of 
parameterization is important. Using the "conventional" parameterization 
by the x coordinate, the largest cross-section appears at x = 3/8, and it 
is equal to 1/2 so that // = 7r 2 /(l/2) 2 « 39.48. Only the first eigenvalue 
Ai is below [i so that the exponential estimate is not formally applicable to 
other eigenfunctions. At the same time, it is clear that the circular branch 
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Figure 9. (Color online) The squared Z/2-norm, J n (xo), of 
two eigenfunctions with n = 1,14 (symbols) for the tilted 
branch on Fig. 2g. The formal estimate (31) with fi = 
7T 2 / (1/4) 2 is plotted by solid line for n = 1. The improved es- 
timate with [i = 7r 2 /(l/4/v / 2)) 2 is shown by dashed (n = 1) 
and dash-dotted (n = 14) lines. 

of "true" width b can be naturally parameterized by the angle (f> or, equiva- 
lently, by an arc, as illustrated on Fig. 2h. However, there is an ambiguity 
in the choice between arcs of various radii (e.g., inner, outer or middle arcs). 
Although the length of all these arcs is proportional to the angle 4>, the pro- 
portionality coefficient enters in the decay rate. For the results plotted on 
Fig. 10, the inner arc of radius r = 3/8 was used. The squared L2-norm was 
plotted as a function of the curvilinear coordinate x' = r(f>, with (ft varying 
between and ir/2. For such a curvilinear parameterization, the width of 
the branch is constant, b = 1/4, so that there are 9 eigenvalues A n below 
/i = 7r 2 /(l/4) 2 . For instance, the exponential decay of the 9th eigenfunc- 
tion is confirmed on Fig. 10. For circular branches, rigorous estimates can 
be derived by writing explicit series representations of eigenfunctions in a 
circular annulus and reformulating the analysis similar to that of Sec. 2 
(this analysis is beyond the scope of the paper). In general, an extension 
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Figure 10. The squared L2-norm, J n (x' ), of four eigenfunc- 
tions with n = 1,9,10,11 (symbols) for the circular branch 
on Fig. 2h. The estimate (31) with [i = 7r 2 /(l/4) 2 is plotted 
by solid (n = 1) and dashed (n = 9) lines. The curvilinear 
coordinate x' Q = rep with r = 3/8 and the angle <p varying 
from to 7r/2, was used. 

of the derivation in Sec. 3 to curvilinear parameterizations is an interesting 
perspective. 



Branch with a small broadening 



A small broadening in the middle is another interesting perturbation of 
the rectangular branch (Fig. 2i). Although this perturbation is small, the 
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threshold value is reduced to \i = ir /(1/2) 2 « 39.48, instead of the value 
fj,' = 7r 2 /(l/4) 2 fa 157.91 for the rectangular branch. As a consequence, the 
sufficient condition X n < fi is satisfied only for n = 1, while for other n, the 
exponential estimate (13) cannot be applied. One may thus wonder how do 
these eigenfunctions behave in the branch? 

We checked that the 8 first eigenfunctions exponentially decay along the 
branch, similarly to the rectangular case. In turn, the 9th, 10th and some 
other eigenfunctions do not exhibit this behavior. Figure 11 illustrates this 
result for three eigenfunctions with n = 1,8,9. Since the rigorous estimate 
is not applicable, we plot the function J„(0) exp(— lyj \i' — A n xo) with // = 
7r 2 /(l/4) 2 as for the rectangular branch. One can see that this function 
correctly captures the behavior of J n (xo) but fails to be its upper bound 
(there are some regions in which J n (xo) exceeds this function). 

In summary, in spite of the fact that the condition A < [i is not sat- 
isfied, the presence of a small broadening does not significantly affect the 
exponential decay of the first eigenfunctions but the upper bound is not 
valid. 

Bifurcating branch 

After describing the eigenfunctions in a single branch, one may wonder 
about their properties in a general pore network which consists of basic 
domains ("pores") connected by branches. In spite of numerous potential 
applications for diffusive transport in porous media, very little is known 
about this challenging problem. A dumbbell domain with a thin "handle" is 
probably the most studied case [26-29]. In particular, the method of match- 
ing asymptotic series yields the asymptotic estimates when the diameter of 
the channel is considered small parameter. 

In order to illustrate the related difficulties, we consider a rectangular 
branch which bifurcates into two rectangular branches (Fig. 2j). For this 
example, the largest cross-section length is 1 so that fi = tt 2 ~ 9.87 and 
the exponential estimate is not applicable. Figure 12 shows J n (xo) for 
three eigenfunctions with n = 1,7,8 and the hypothetical estimate with 
fj,' = 7r 2 /(l/4) 2 as for the rectangular branch. The first eigenfunction is well 
estimated by the exponential function Ji(0) exp(— — XiXq) for practi- 
cally the whole length of the branch (xq < 1), with a small deviation at the 
bifurcation region (xq > 1). Similar behavior is observed for other eigen- 
functions up to n = 7. The larger the index n, the earlier the deviation 
from the exponential estimate appears (e.g., for n = 7, the estimate works 
only for x$ < 0.3). In turn, the 8th eigenfunction has no exponential decay 
in the branch as it is localized in the bifurcation region. At the same time, 
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Figure 11. (Color online) The squared L2-norm, J n (xo), 
of three eigenfunctions with n = 1,8,9 (symbols) for the 
branch with a broadening on Fig. 2i. The estimate (31) 
with fi = 7r 2 /(l/2) 2 for n = 1 is shown by dotted line. The 
hypothetical estimate (31) with jjl = 7r / (1/4) 2 is plotted by 
solid (n = 1), dashed (n = 8) and dash-dotted (n = 9) lines. 
Vertical dotted lines indicate the position of the broadening. 



the corresponding eigenvalues are close: A7 ~ 125.20 and As ~ 127.67. This 
illustrates how significant may be the difference in properties of two consec- 
utive eigenfunctions. We also emphasize on the difficulty of distinguishing 
these cases in general. 
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Figure 12. (Color online) The squared L2-norm, J n (xo), 
of three eigenfunctions with n = 1,7,8 (symbols) for the 
bifurcating branch on Fig. 2j. The hypothetical estimate 
(31) with \i = 7r 2 /(l/4) 2 is plotted by solid (n = 1), dashed 
(n = 7) and dash-dotted (n = 8) lines. The vertical dotted 
line indicates the position of bifurcation. 

Surprisingly, an exponential decay can be recovered for thinner bifurca- 
tions of arbitrary length. Although this example may look specific, the 
arguments provided below are of general interest. In order to prove the ex- 
ponential decay of the squared norm J(xq) in Eq. (29), we will establish the 
inequalities similar to those for I(xq) of Sec. 3. Replacing the norm I(xo) 
by the weaker norm J{xq) allows one to produce more accurate estimates. 
Although we bear in mind the example shown on Fig. 2j , the arguments are 
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rather general and applicable in higher dimensions. By definition, one has 
J'{x Q ) = -I(x ) = - J u 2 (x ,y)dy <0 



tt(x ) 

and 

J"(xq) = — 2 J u—udy = 2 j [(V«, Vu) + uAji] dxdy, 

Q(xo) Q(xo) 

where the Green's formula, the identity d/dx = —d/dn on £l(xo) and Dirich- 
let boundary condition on dQ(xo)\Q,(xo) were used. The first integral can 
be estimated by splitting Q(xq) into two subdomains: Q\(xq) (a part of the 
rectangular branch) and Qi (the bifurcation): 

J (Vu, Vu)dxdy > [i J u 2 dxdy + J (Vu, Vu)dxdy, 

Q(xo) Qi(xo) Q2 

where the Friedrichs-Poincare inequality was used for Qi(xq), with the small- 
est fj, over Qi(xq) (for the rectangular branch on Fig. 2j, \x = -K 2 /b 2 ). 
The separate analysis of two subdomains allowed us to exclude large cross- 
sections of the bifurcation Q2 from the computation of /i. 

In order to estimate the second integral over Q2, we note that 

. (Vv,Vv) 
vi = mi — r— 

vGH 1 (Q 2 ), v\ r =0, v^O (v,V) 

is the smallest eigenvalue of the Laplace operator in the rectangle Q2 with 
Dirichlet boundary condition on the top, bottom and right segments (de- 
noted by r) and Neumann boundary condition on the left segment. Since 
the eigenfunction u also belongs to H 1 (Q2), one has 

(Vu,Vn) L2(Q2 ) > vi(u,u) L2(Q2) , 

from which 

J" (a? ) > 2(/i - A) J u 2 dxdy + 2(vi - A) J u 2 dxdy. 

Qi(xo) Q2 

If v\ > fx, we finally obtain 

J"{xq) > 2(ii - A) J u 2 dxdy = 2(ii- \)J(x ). 

Q{xo) 

As shown in Sec. 3, this inequality implies an exponential decay of J(xq). 
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If Q2 is a rectangle of height h and width w with Dirichlet boundary con- 
dition on the top, bottom and right segments and Neumann boundary con- 
dition on the left segment, the smallest eigenvalue is v\ = it 2 / h 2 + tt 2 / (2w) 2 . 
If the branch Q\ is also a rectangle of width 6 (as shown on Fig. 2j), one 
has fi = TT 2 /b 2 . The condition v\ > fi reads as (b/h) 2 + (b/(2w)) 2 > 1. 
For instance, this condition is satisfied for any h if the bifurcation width 
w is smaller than 6/2. In particular, this explains that a small broad- 
ening of width 6/2 shown on Fig. 2i does not degrade the exponential 
decay of J n (xo) for xo < 0.4375. In turn, the bifurcation shown on Fig. 
2j with h = 1 and w = b = 1/4 does not fulfill the above condition as 
(b/h) 2 + (b/(2w)) 2 = 1/16 + 1/4 < 1. 

4.2 Infinite branches 

In the derivation of Sec. 3, there was no condition on the length a of the 
branch. Even for an infinite branch, the exponential decay is valid once the 
condition (11) is satisfied. However, when the branch is infinite, the Laplace 
operator eigenspectrum is not necessarily discrete so that L2-normalized 
eigenfunctions may not exist. A simple counter-example is a semi-infinite 
strip D = [0, 00) x [0, 7r] for which functions sinh(\/^ 2 — Ax) sin(ny) satisfy 
the eigenvalue problem (8) but their Li2(D) norms are infinite. 

As shown by Rellich, if the first eigenvalue fi\(x) in the cross-section 
Q(x) of an infinite branch Q goes to infinity as x — > 00, then there ex- 
ist infinitely many /^-normalized eigenfunctions [30]. In two dimensions, 
Hi(x) = n 2 /£(x) 2 is related to the length l(x) of the largest interval in the 
cross-section £l(x). The condition /^i(x) — > 00 requires thus i(x) — > 0. At 
the same time, the number of intervals in O(x) can increase with x so that 
the "total width" of the branch may arbitrarily increase. 

For infinite decreasing branches, eigenfunctions can be shown to decay 
faster than an exponential with any decay rate. In fact, although the thresh- 
old [i was defined in (10) as the smallest \i\{x) over all cross-sections of the 
branch, the separation into a basic domain and a branch was somewhat 
arbitrary. If the separation occurs at xq (instead of 0), one gets 

I(x) < I(x ) exp[-2 y 7 /x(x ) - A (x - x )] x > x , 
where the new threshold h(xq) = inf [ii(x) increases with xq, while the 

xo<x 

prefactor I(xq) also decays exponentially with xq according to Eq. (13). 
Since the above estimate is applicable for any xq, one can take xq to be a 
slowly increasing function of x [its choice depends on /x(xq)] that would result 
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in a faster-than-exponential decay of I(x). This result is a consequence of 
the condition fii(x) — > oo as x goes to infinity. 

4.3 Three-dimensional domains 

The exponential estimate (13) becomes still more interesting in three (and 
higher) dimensions. While any cross-section Q(x) was a union of intervals 
in two dimensions, the shape of cross-sections in three dimensions can vary 
significantly (e.g., see Fig. lc). Whatever the shape of the branch is, the only 
relevant information for the exponential decay is the smallest eigenvalue \i in 
all cross-sections Q(x). For instance, for a rectangular profile of the branch, 

2 2 

tt(x) = [0, b(x)] x [0,c(x)], the first eigenvalue fii(x) = + can 
remain bounded from below by some \x even if one of the sides b{x) or c(x) 
grows to infinity. This means that eigenfunctions may exponentially decay 
even in infinitely growing branches. 

4.4 Neumann boundary condition 

The theoretical derivation in Sec. 3 essentially relies on the Dirichlet bound- 
ary condition on the branch boundary: u\qq = 0. This condition can be 
interpreted, e.g., as a rigid fixation of a vibrating membrane at the bound- 
ary, or as a perfect absorption of diffusing particles at the boundary. The 
opposite case of free vibrations of the membrane or a perfect reflection of 
the particles is described by Neumann boundary condition, du/dn\QQ = 0. 
Although the eigenvalue problem may look similar, the behavior of eigen- 
functions is different. In particular, an extension of the results of Sec. 3 fails 
even in the simplest case of a rectangular branch, as illustrated on Fig. 13. 
Although the 8 first eigenvalues A n are below fi = 7r 2 /(l/4) 2 , only some of 
them decay exponentially (e.g., with n = 4). This decay seems to be related 
to the reflection symmetry of the domain. Since the eigenfunctions with 
n = 4, 10 are anti-symmetric, they are along the horizontal symmetry line. 
One can therefore split the domain into two symmetric subdomains and im- 
pose the Dirichlet boundary condition on the splitting line. Although an 
exponential decay may be expected for this new problem, its mathematical 
justification is beyond the scope of the paper. 

On the other hand, the derivation of the exponential estimate in Sec. 3 
does not use the boundary condition imposed on the basic domain V. As 
a consequence, the estimate is applicable for arbitrary boundary condition 
on V which guarantees the Laplace operator in the whole domain D to be 
self-adjoint. 
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n = 2 n = 3 




Figure 13. Six eigenfunctions with n = 2,3, 4, 5, 9, 10 for 
the rectangular branch (Fig. 2a) with Neumann boundary 
condition (the fundamental eigenfunction with n = 1 is con- 
stant and not shown). 

4.5 Expelling from the branch 

It is important to stress that the "smallness" of an eigenfunction in the 
branch and its exponential decay are different notions which should not to 
be confused. In fact, the eigenfunction can be small in the branch either 
due to a rapid exponential decay, or because of the small constant I n (0) or 
J n (0) in front of the estimate. For instance, Fig. 10a shows the behavior of 
^10(^0) without an exponential decay, but the eigenfunction is nevertheless 
small (Fig. 10b). Another example on Fig. 4a with n = 9 illustrates the 
opposite situation: the eigenfunction decays exponentially along the branch 
but it does not look small. 
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Conclusion 

We have studied the behavior of the Laplace operator eigenfunctions in a 
large class of domains composed of a basic domain of arbitrary shape and a 
branch Q which can be parameterized by a variable profile Q(x). We have 
rigorously proved that each eigenfunction whose eigenvalue A is smaller than 
the threshold fi = inf{/Ui(x)}, exponentially decays inside the branch, where 
H\{x) is the first eigenvalue of the Laplace operator in the cross-section Vt(x). 
In general, the decay rate was shown to be at least \[2\J \x — A. For non- 
increasing branches, the decay rate 2\/ [i — A was derived and shown to be 
sharp for an appropriate parameterization of the branch. The exponential 
estimate is applicable in any dimension and for finite and infinite branches. 
In the latter case, the condition \x\{x) — > oo as x — > oo is imposed to ensure 
the existence of L2-normalized eigenfunctions. Since the derivation did not 
involve any information about the basic domain V, the exponential estimate 
is applicable for arbitrary V with any boundary condition on dV for which 
the Laplace operator in D is still self-adjoint. In turn, the Dirichlet boundary 
condition on the branch boundary was essential. Note that the mathematical 
methods of the paper can be adapted for studying eigenfunctions for various 
spectral problems or other types of domains. 

The numerical simulations have been used to illustrate and extend the 
theoretical results. It was shown that the sufficient condition A < /i is 
not necessary, i.e., the eigenfunctions may exponentially decay even if A > 
ji. However, in this case, the decay rate and the range of its applicability 
strongly depend on the specific shape of the branch. For all numerical 
examples, the sharp decay rate 2y/ fj, — A was correct, even if the condition 
(15) for non-increasing branches was not satisfied. In future, it is tempting 
either to relax this condition, or to find counter-examples, for which the 
sharp decay is not applicable. 
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From the inequality sinhx < coshx, Eq. (5) is bounded as 

b °° T vr 1 /" 

\\ Vu \\l 2 (Q(x )) < 2 12 C n (^ n ) 2 + ^n COsll 2 (~f n (a-x))dx. 

The last integral is estimated as 

a a 

j cosh 2 ( 7n (a - x))dx < j e 2 ^ a ~ x ^dx = 



3 27„a p 2-y n a p 2-y n a 



271x0 



27n 27„ 27„ 

where we used the inequality coshx < e x for x > and the fact that 
7 n = y / TT 2 n 2 /b 2 — A increases with n. We have then 



h 00 
l|V«llL ( g ( , ))<^- 2 ^°E C n 



2 



2 2 
6 2 2 7n + 2 



e 27na . 



Writing two inequalities: 



7n = ^Jir 2 n 2 /b 2 - A < ^n, 



n 2 n 



7n (ir/b)y/l- Xb 2 /(TT 2 n 2 ) 



< Cm, 



where C\ is a constant, one gets an upper bound in the order of n for the 
expression in large brackets. Finally, we have an estimate for e 27na as 

sinh 2 (2 7n a) = ^(1 - e" 2 ^) 2 > ^(1 - e" 2 ^) 2 , 
from which 

e 2 ^ a <C 2 sinh 2 (2 7n a), 

with a constant C2 = 4/(1 — e _27ia ) 2 . Bringing together these inequalities, 
we get the estimate (6). 



Trace theorem 

The trace theorem implies [31] that the series 

,/hY, 

'L 2 (0,&)> 



f( x ) = ^2n{u(x,y),sm(irny/b)) 2 } 



n=l 
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1/2 

which is equivalent to the squared norm of u(x, y) in the Sobolev space HJ ^ , 
may be estimated from above by the norm ||Vu||^ 2 ^. For completeness, 
we provide the proof for our special case. 

For a fixed x, we denote 

X n {x) = (u(x,y),sm(irny/b)) L2{ob) 

the Fourier coefficients of the function u(x,y): 



2 x \ 

u(x, y) = - X n {x) sin(vrny/6). 



n=l 



On one hand, starting from X n (a) = 0, one gets 





a 




a 


X 2 n (x) = 


j (x 2 Jd Xl 


= 2 


J X n X' n d Xl 




X 




X 



while the Cauchy inequality implies 
X n X' n dx\ 



< 2||X n || i2 ( 0)a )||X n || i2 ( 0jO ). 



The inequality 2af3 < a 2 + /3 2 yields 

2(7rn/6)||X n || La ( 0iO )||XllL a (0,o) < ll X nlll 2 (0,a) + ( W 6 ) 2 | l^nl ll 2 (0,a) > 
from which 

(™/b)xfa) < nx;ni 2(0iO) + (w&) 2 n^niiL (0 , a) . 

On the other hand, we write explicitly the energetic norm of u: 



|Vu| 



L 2 (Q) 



2 00 / \ 

^E(lKllL(0,a) + (W6) 2 |l^n|li 2(0 , a) j: 
n=l ^ ' 



from which 



n=l 



f{x) = J2nX*(x) < — ||Vu| 



2 < _||V7/H 2 

L 2 (Q) S 2n \\VU\\L 2 (D)- 



Since the coefficients X n {x} and. c n are related as 



X n {x) = - c n sinh(7„(a - x)), 
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the substitution of x = into this equation yields 
(32) 

4 4 2A 2A 

^n C 2sinh 2 ( 7n a) = p^^n(0) 2 = p/(0) < — \\u\\l 2{D) = — |M|£ a(r)) . 

n=l n=l 



Appendix B Several classical results 

For completeness, we recall the derivation of several classical results [31, 32] 
which are well known for spectral analysts but may be unfamiliar for other 
readers. 



B.l Rayleigh's principle 

Let us start with the first eigenvalue Ai of the problem (1) which can be 
found as 

(33) A, = inf ^ VV ^\ 

O 

where H 1 = {v G L 2 {D), dv/dxi £ L 2 (D), i = l,...,n + 1, v\ dD = 0}. 
Denoting the first eigenfunction in Eq. (28), one takes 

v= Uu {x,y)eV, 
[0, (x,y)$V, 

as a trial function in Eq. (33) to obtain 

(V<£i, V0i) i2 (v) 

m,(pi)L 2 (v) 

i.e., the first eigenvalue Ai in the whole domain D is always smaller than 
the first eigenvalue «i in its subdomain V. More generally, if there are n 
eigenvalues «i < • • • < n n < /i then there exist n eigenvalues Ai < ■ ■ ■ < 
A„ < fi. 

Note that the Friedrichs-Poincare inequality (21) follows from (33). 



B.2 Rellich's identity 



Let u be an eigenfunction which satisfies the equation 

Au + Xu = (x,y)eD, u\ 9D = 0. 
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We multiply this equation by |^ and integrate over the domain Q(xq) defined 
by Eq. (23): 

/3u f dix 

—An dxdy + A / u— dxdy = 0. 

Q(xo) Q(xq) 

The second integral can be transformed as 



A / u d ^dxdy=\ J [^dxdy 

Q{x ) Q(x ) 
( 35 ) X r X r 

/ u 2 (x ,y)(e x ,n)dS = -- / u 2 (x ,y)dy, 



2 

8Q(xo) n(aio) 

where n = n(S) is the unit normal vector at S £ dQ(xo) and the boundary 
condition uqd = was used on r(xo) = dQ(xo)\£l(xo). 

Using the Green's formula, the first integral in Eq. (34) can be trans- 
formed to 

,36) f^y.f^JS- f(vu,V^y. 

Q(x ) 8Q(x ) Q(x ) 

The first integral over dQ(xo) can be split in two terms: 

dudu^_ I" d 11 ® 11 ^^ f ( ' 9u\ ^ 
dx dn J dx dn J \dx J ' 

OQ(x ) T(x ) U{x ) 

where d/dn = (n, V) is the normal derivative pointing outwards the domain, 
and the sign minus appears because du/dn = —du/dx at Q,{xq). 

The second integral in Eq. (36) is 

' du\ 1 f d 

Vu, V— dxdy = - / —(Vu,Vu)dxdy 
K ox J 2 J ox 

Q(xo) Q(xo) 

= - J {Vu,Vu)(e x ,n)dS + - J (Vu,Vu)dy. 

r> ) n(x ) 

Taking into account the Dirichlet boundary condition u|r( x - ) = 0, one has 

it-, m ® u \ 
(V«)|r(x ) =n^|r(x ), 

du \du 
fahzo) = ( e x.V«) = {e x ,n) — \ r(xo) . 



n 
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Combining these relations, one gets 




Q(xo) 



r(x ) 



2 



1 




(V±u, V±u)dy, 



from which and Eqs. (34, 35) the Rellich's identity (22) follows. 
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